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Average phase factor in the PNJL model 
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The average phase factor {e' 2 ' e } of the QCD determinant is evaluated at finite quark chemical potential (/x q ) 
with the two-flavor version of the Polyakov-loop extended Nambu-Jona-Lasinio (PNJL) model with the scalar- 
type eight-quark interaction. For /i q larger than half the pion mass m w at vacuum, (e 2 ' e ) is finite only when the 
Polyakov loop is larger than ~ 0.5, indicating that lattice QCD is feasible only in the deconfinement phase. A 
critical endpoint (CEP) lies in the region of (e 2! ) = 0. The scalar-type eight-quark interaction makes it shorter 
a relative distance of the CEP to the boundary of the region. For /i q < m«/2, the PNJL model with dynamical 
mesonic fluctuations can reproduce lattice QCD data below the critical temperature. 

PACS numbers: ll.30.Rd, 12.40.-y 



I. INTRODUCTION 

The thermodynamics of quantum chromodynamics (QCD) 
is well defined, since QCD is renormalizable and parameter 
free. Nevertheless, the thermodynamics is not understood 
at lower temperature (T) because of its nonperturbative na- 
ture. The thermodynamics is closely related to not only natu- 
ral phenomena such as compact stars and the early universe 
but also laboratory experiments such as relativistic heavy- 
ion collisions. Lattice QCD (LQCD) is a first-principle cal- 
culation, but it has the well-known sign problem when the 
quark-number chemical potential (/i q ) is real; for example, 
see Ref. QJJ] . So far, several approaches have been proposed 
to circumvent the difficulty; for example, the reweighting 
method [2], the Taylor expansion method yfl and the analytic 
continuation from imaginary /j, q to real fi q Jj-0]. However, 
these are still far from perfection particularly at /J, q /T > 1. 

The success of the approaches is linked to how difficult the 
sign problem is. As a good index of the difficulty, one can 
consider the average of the phase factor 



j.ie = det(£> + /z q 7o + to) 
det(D + ,u q 7o + to)* 



(1) 



of the Fermion determinant. If the average of the phase factor 
is much smaller than 1, this means that there are severe can- 
cellations in the path integral of the QCD partition function. 
In this situation, LQCD simulations are not feasible. 

The average is obtained by taking the expectation value of 
the phase factor in the phase-quenched theory in which the 
Fermion determinant is replaced by the absolute value. In the 
two-flavor case, the average is 



(e 2W ) = 



n+i 



Zi 



(2) 



+i* 
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where Zi+i stands for the partition function of the ordinary 
two-flavor theory and Zi+i* represents that of the two-flavor 
phase-quenched theory in which one of two flavors is changed 
into a conjugate flavor. For comparison of the 1+1* system 
with the 1+1 system, let us introduce the modified isospin 
chemical potential /ij related to the isospin chemical poten- 
tial /ij so as /it = ^i so /2. When the 1+1 system has a value of 
Hq, the 1+1* system possesses the same value of /Ut. 

It is not easy to calculate the average phase factor with 
LQCD even for small (J>q/T. Actually, several LQCD results 
on the average phase factor are spotted; see Ref. [11] and ref- 
erences therein. It is then important to make a systematic anal- 
ysis on the phase factor by using effective theories. This was 
done by the chiral perturbation theory [i§, |9fl. The result is con- 
sistent with LQCD one [11] when T is lower than the critical 
one T c . However, the theory is not valid for T > T c . 

Recently, the average phase factor at T > T c was calcu- 
lated by the random matrix theory [ 10] and the Nambu-Jona- 
Lasinio (NJL) model [12]. When the saddle-point approxima- 
tion is applied to the path integrals in the partition functions 
Zi+i and Zi+i* , the average phase factor can be described by 



„2iB\ 



y/detH: 



1+1* 



>/det Hi 



-j9V(Oi 



-o. 







(3) 



+ i 



where j3 — 1/T, fl is the thermodynamic potential at 
mean field level and H is the Hessian matrix showing 
static fluctuations (SF) at the saddle point. The average 
phase factor thus obtained is dominated not by the ex- 
p onential fac tor e ~^ v( - ai +1 ~ ni + 1 *'> but by the SF factor 
y'det i?i+i* / y'det -ffi+i IU2I1 . because J7i+i = i?i+i» in 
the normal phase with no pion condensate and the SF factor 
is zero in the pion condensate phase; this will be explained 
explicitly in subsection MB I Thus, the average phase factor 
should be calculated with the mean field (MF) approximation 
plus SF corrections. This framework is referred to as MF+SF 
in the present paper. 

The NJL model describes the chiral phase transition IU3I - 
Ual . but not the deconfinement transition. The Polyakov- 
loop extended Nambu-Jona-Lasinio (PNJL) model [ 19-38] 
was constructed to treat both the transitions simultaneously. 
Very recently, the PNJL model was shown to be successful 



in reproducing LQCD data for imaginary quark chemical po- 
tential jH,|35[], imaginary quark and isospin chemical poten- 
tials Oal and real isospin chemical potential ID8I1 . Thus, the 
PNJL model is one of the most reliable models. 

In this paper, we evaluate the average phase factor by the 
PNJL model in the MF+SF framework and investigate a re- 
lation between the Polyakov-loop and the average phase fac- 
tor. For thermal systems at /ii = /i q = and at /ii > 
and fi q = 0, the scalar-type eight-quark interaction is essen- 
tial for the PNJL model to reproduce LQCD data fH [38^1 . 
We then analyze an effect of the eight-quark interaction on 
the average phase factor, too. Finally, we consider dynamical 
mesonic fluctuations (DF), instead of static mesonic fluctua- 
tions, to compare the PNJL results with LQCD data. 

This paper is organized as follows. In Sec. II, we recapit- 
ulate the PNJL model and a way of calculating the average 
phase factor. The numerical results are shown in Sec. III. Ef- 
fects of the Polyakov loop, the eight-quark interaction, static 
fluctuations and dynamical mesonic fluctuations are investi- 
gated. Section IV is devoted to summary. 



II. FORMALISM 

In this section, we first explain the PNJL model in the MF 
level for the case of finite fj, q and ^i and treat SF in order to 
evaluate the average phase factor. 



For later convenience, we use fi\ instead of ^; so and call it 
the isospin chemical potential simply. The 1+1 and the 1+1* 
theory correspond to taking (^ q ,^i) = (/i q ,0) and (0,/j, q ), 
respectively, in C of (0). In the limit of too = fij = 0, the 
PNJL Lagrangianhas the SU L (2) x SU R (2) x U v (l) x SU C (3) 
symmetry. For too ^ and /ii ^ 0, it is reduced to Ui 3 (1) x 
J7 V (1) x SU C (3). 

The Polyakov-loop operator <P and its Hermitian conjugate 
?2>t are defined as 



<Z> 



l T ^ 



(p* 



i Tri » 



with 



I/(x) = V exp 



i I cItA4(x,t) 
o 



(7) 



(8) 



where V is the path ordering and A4 = iA°. In the Polyakov 
gauge, L can be written in a diagonal form in color space [20]: 

£ _ e «/3(<feA3+08A 8 ) /g\ 

In the MF level, ^3 and <j>s are treated as classical vari- 
ables JH. 

The spontaneous breakings of the chiral and the Uj 3 (l) 
symmetry are described, respectively, by the chiral conden- 
sate a = (qq) and the charged pion condensate [28] 



A. PNJL model 

The Lagrangian of the two-flavor PNJL model in Euclidean 

spacetime is 

£ = q{lvD v - 74A + m )q + G s [{qq) 2 + (qi^ 5 fq) 2 ] 

+ G sS [(qq) 2 + [qi^rq) 2 } 2 -U^[A},$[A)* ,T), 

(4) 

where D v = d v + iA v and A v = 5^gA° a ^- with the gauge 
field A v a , the Gell-Mann matrix A a and the gauge coupling 
g. The coupling constant G s (G s s) represents a strength 
of the scalar-type four-quark (eight-quark) interaction. The 
Polyakov potential U of (fT8l is a function of the Polyakov 
loop <P and its Hermitian conjugate <P*. 

The chemical potential matrix (i, is defined by ft = 
diag(^ u , Hd) with the w-quark (<i-quark) number chemical po- 
tential /i„ ip.d), while mo = diag(mo,TOo). The chemical 
potential matrix is described by /i q and /ij as 



IX = n q T + nm, 



(5) 



where tq is the unit matrix and n (i = 1, 2, 3) are the Pauli 
matrices in flavor space. The /i q and [i\ are related to the 
baryon and isospin chemical potentials, /jb and /ij SO , coupled 
respectively to the baryon charge B and to the isospin charge 
7/ 3 as 



Mq 



Mu + Hd 



MB 

T' w = 



fJ>u ~ Vd 



~2 



= SF- (6) 



tt ± = ^e ±M = (qij 5 T±q). 



(10) 



Since the phase a represents the direction of the U\ 3 (1) sym- 
metry breaking, we take a — for convenience. The pion 
condensate is then expressed by 



7T = (qi'ysriq). 



(ID 



Making the MF approximation 1118112811 . one can obtain the 
MF Lagrangian as 

£mf = q{li>D u - j4fj, + Mt + Nij 5 Tx)q 

+G s [<T 2 + 7r 2 }+3G s8 (<T 2 +ir 2 ) 2 +U (12) 

with 

M = too - 2[G S + 2G s8 (a 2 + n 2 )}a, 

N = -2[G s + 2G s8 (<T 2 +ir 2 )}ir. (13) 

The quark propagator S in the MF level, that has off-diagonal 
elements in flavor space, is obtained by 



s-'ip) 



1»P V - M«74 + M 



-i l5 N 
luP u - Md74 + M 



(14) 



where P v = p v — A v . Performing the path integral in the 
PNJL partition function 



Z 



MF 



DqDq exp 



<TxC 



MF 



(15) 



we can get the thermodynamic potential (per unit volume), 



B. Average phase factor 



n MF = -Tln(Z MF )/V = -2j2 



i=± 



d 3 p 

(2tt) 3 



3^(p) 



with 



~ In [1 + 3(* + $T^"W)e-^"fp) + e~ 3 ^r(p)] 
- In [1 + 3(<T + <?> e -^ l + (p))e-^, + (p) + e -3^+(p)]- 
G s [a 2 + t: 2 } + 3G s8 {a 2 + t: 2 ) 2 + U 

^(p)=^±(p)±/i q , 



(16) 



(17) 



where E± (p) 



V(^(p)±w) 2 +^ 2 and £(p) = 



\/p 2 + M 2 . The mometum integral in (fTBl is regularized by 
the three-dimensional cutoff A 
We use W of Ref. |J25ll : 



U = T 4 



a(T) 



<P*$ 



b(T) ln(l - 6<2><£* + 4(<£ 3 + <P* 

'T \ , /T \ 2 



3(<£<T) 2 ) 



(18) 



Tn\3 



a (r) = ao+a 1 (^)+a 2 (^) , 6(T) = 6 3 (y 



(19) 



The parameters of W are adjusted to LQCD data in the heavy- 
quark (pure-gauge) limit [39, 40]; the resultant parameter set 
is shown in Table I. In the limit, the Polyakov potential yields 
a first-order deconfinement phase transition at T — Tq. Since 
the first-order transition takes place at T = 270 MeV in 
LQCD, T is often set to 270 MeV in the PNJL calculation. 
In the light-quark case, however, the PNJL calculation with 
this value of To yields a larger value of T c than the full LQCD 
prediction T c = 173 MeV MlrSl- Therefore, we rescale T 
to 200 MeV to reproduce T c = 173 MeV iQ. 



«o 



«i 



(12 



h 



3.51 



-2.47 15.2 



-1.75 



TABLE I: Summary of the parameter set in the Polyakov-potential 
sector determined in Ref. 1 25]. All parameters are dimensionless. 



In the NJL sector, two parameter sets are taken; in the 
first set G s8 is finite, while in the second set it is zero. The 
first set is A = 0.6315 GeV, G s = 4.673 [GeV~ 2 ], G s8 = 
452.12 [GeV -8 ] and mo — 5.5 MeV. This set reproduces not 
only the pion decay constant / w = 93.3 MeV and the pion 
mass m-x = 139 MeV at vacuum but also LQCD data [41-43] 
on a and |«P| for thermal systems with no /i q and /ii [35]. The 
second set is A = 0.6315 GeV, G s = 5.498 GeV" 2 , G s8 = 
and mo = 5.5 MeV. This set can reproduce the pion decay 
constant and the pion mass at vacuum, but not LQCD data for 
thermal systems with no /z q and \i\. Thus, the first set with 
finite G s $ is more reliable. 



As mentioned in sectionU we have to consider fluctuations 
to mean fields to evaluate the average phase factor. In this 
subsection, we consider static fluctuations (SF). 

In the path-integral representation of the partition function 
Z or the thermodynamic potential J?, <fi 3 and </) 8 are funda- 
mental fields rather than <P and <P* . This means that we should 
solve the stationary conditions 



df2 

dip 







(20) 



for if = (ct,7t, 03,(/> 8 ) rather than ip = (a, if, <?,<?*) B25I1 . 
although the solutions are not so different between the two 
cases. However, the first case does not guarantee that Q is 
real. Following Ref. [25], we then put 0g = to keep Q real. 
Noting that the first-order derivative with respect to 0§ does 
not vanish, we expand J? up to quadratic terms of fluctuations: 



n = n 



MF 



6f2\ i / s 2 n 

8(fiJ UF l 2 \8ipi6tpjs MJ , 



SifiSifj, (21) 



where ip = J2i ^fi + <£>mf for mean fields ipMF and static 
(constant) fluctuations Sipi. Since first-order terms in Sipi are 
purely imaginary, we can regard an integral over Sipi as a 
Fourier integral. We then obtain 



Z = JJd(<5y>j)exp 



V 

T' 



n 



i 



exp 



T 



with 



and 



1 f 5 2 Q S ~' 



1 2 VO^08/MF VO08/MF 



A-=(^)'||<totff||*, 



(22) 



(23) 



(24) 



where n is the number of fields and H is the Hessian matrix 
defined by 



H 



s 2 n 



SipiSifj 



(25) 



MF 



Obviously, the Hessian matrix describes static fluctuations of 
mean fields <^mf- Therefore, the average phase factor is ob- 
tained by ([3} with Q replaced by Q. The average phase factor 
thus obtained is a function of the ipuF that satisfy the station- 
ary conditions (1201 1, Thus, the average phase factor is calcula- 
ble in the MF+SF framework. 

The static fluctuations are composed of those of a and it and 
of 03 and 0g. In subsection lHICl we keep treating static fluc- 
tuations for 03 and 8 , but consider dynamical fluctuations for 
a and -k, that is, a and 7r mesons. 

It was revealed in Ref. |[12ll with the NJL model that in the 
pion condensate phase, where 7Tmf is finite and then a mass- 
less mode appears, the average phase factor vanishes owing 



to dot H\ + i* = 0. Obviously, this is true also for the PNJL 
model, as shown in ( f25t . In the normal phase with no pion 
condensate, J?i+i and J?i+i* are the same in the MF level, 
so that the average phase factor is determined by only the SF 
factor ^/det i?i+i« /ydet #1+1 • 

LQCD calculation in Ref. IQ] has a lattice size 16 3 x 4. 
Hence, the three-dimensional volume is V — (16a) 3 for a lat- 
tice spacing a and the inverse of temperature is 1/T = 4a. 
The four-dimensional volume is then obtained by V/T = 
64T . This four-dimensional volume is taken also for the 
PNJL model in its calculation of the average phase factor. 



III. NUMERICAL RESULTS 

A. Average phase factor and Polyakov loop 

Solving the stationary conditions ( f20b and inserting the so- 
lutions in (01, one can obtain <P, it and (e 2i8 ) as a function of 
T and /j, q . As a shorthand notation, we use <P for <?i+i and 
it for 7Ti+i*. All the calculations in this subsection are done 
without the scalar-type eight-quark interaction; roles of the 
eight-quark interaction will be discussed in subsection lHIBI 
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Fig. 1 : <P dependence of the scaled pion condensate n/ao and the 
average phase factor (e 2ie ) at /i q = 100 MeV. Here, ao is a chiral 
condensate at T = /x q — 0. The former is plotted by a dashed curve, 
while the latter is by a solid curve. 

In Fig. |U we plot (<2>(T),7r(T)) by a dashed line and 
(#(T), (e 2ie )(T)) by a solid line, varying T with //„. fixed 
at 100 MeV. Since <P is an increasing function of T 1 35], an 
increase of <P means that of T. The pion condensate it de- 
creases as <P increases and finally vanishes at a critical value 
<P C of $. Below <P C , the average phase factor is always zero, 
while it is finite. Above <P C , inversely, the average phase fac- 
tor is finite, while it is always zero. Thus, there is a negative 
correlation between the average phase factor and the pion con- 
densate. This property is also seen in the NJL model [12]. In 
contrast, there exists a positive correlation between the aver- 
age phase factor and the Polyakov loop: the average phase 
factor is zero at small <P such as <P < <P C , but at large <P such 
as <P > <P C the average phase factor is finite and an increasing 
function of <P. In the <P = 1 limit, the average phase factor 



tends to 1 . This implies that the NJL model overestimates the 
average phase factor compared with the PNJL model, since 
<P = 1 in the NJL model. This is true, as shown below in 
Fig.0 




KjtGeV] 

Fig. 2: Boundary of the (e 2 * } = region in /i q -T plane. A 
solid (dashed) line represents the boundary in the PNJL (NJL) model. 
Symbols x and + stand for critical endpoints of the PNJL and NJL 
models, respectively. 

Figure|2]shows a boundary of the region where (e 2ie ) = 0. 
The solid (dashed) curve is a result of the PNJL (NJL) model. 
The region is wider in the PNJL model compared with the 
NJL model. Thus, the fact that $ < 1 in the PNJL model 
makes the region wider and the average phase factor smaller 
outside the region. The critical endpoint (CEP) of the first- 
order chiral phase transition is also plotted by a cross (x) for 
the PNJL model and a plus (+) for the NJL model. For both 
the models, thus, the CEP and the first-order phase transition 



are in the 



JliBs 



region. This implies that the location 



of CEP cannot be determined by LQCD directly. The CEP is 
located at higher T and lower /i q in the PNJL model compared 
with the NJL model. However, a relative distance of CEP 
to the boundary of the (e 2,e ) = region is almost the same 
between the two models. 



B. Scalar-type eight-quark inteaction 



The scalar-type eight-quark interaction 

GssKqqf + (^7 5 rg) 2 ] 2 



(26) 



is inevitable to reproduce LQCD data at finite real isospin 
chemical potential IBal and zero and finite imaginary quark 
chemical potential 113511 - Furthermore, the interaction affects a 
location of CEP for real /x q . Therefore, it is important to see 
how the scalar-type eight-quark interaction affects the average 
phase factor. 

Figure [3] shows /i q dependence of the average phase factor 
at T = 0.9T C , T c and 1.1T C . The solid and dashed lines cor- 
respond to the PNJL result with and without the eight-quark 
interaction. Below T c such as T = 0.9T C , the scalar-type 
eight-quark interaction does not affect the average phase fac- 
tor, and at T = T c the effect becomes appreciable. Above 



T c such as T = 1.1T C , it enhances the average phase factor 
sizably. 



proportional to the relative distance. In this sense, the fact that 
the relative distance is small is important. 
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Fig. 3: Effect of scalar-type eight-quark interaction on average phase 
factor. Solid and dashed lines stand for the PNJL result with and 
without the eight-quark interaction, respectively. 

Figure 2]presents a boundary of the (e 2%e ) = region and 
locations of CEP calculated by the PNJL model without and 
with the scalar-type eight-quark interaction. The eight-quark 
interaction makes the region smaller. Meanwhile, it shifts the 
CEP to higher T and lower /i q . Thus, the relative distance 
of CEP to the boundary becomes much smaller by the scalar- 
type eight-quark interaction, although CEP itself lies in the 
(e 2t0 ) = region even after the scalar-type eight-quark in- 
teraction is taken into account. If more accurate LQCD data 
becomes available in future outside the region (e 2lB ) = 0, the 
PNJL model that reproduces the data can predict a location 
of CEP in principle. The reliability of the prediction may be 
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Fig. 4: Effect of the scalar-type 8-quark interaction on locations of 
CEP and the boundary of the (e 2 ' e ) — region. Solid and dashed 
lines stand for the boundary calculated by the PNJL model with and 
without the eight-quark interaction, respectively. Symbols + and x 
mean the CEP calculated by the PNJL model with and without the 
eight-quark interaction, respectively. 

Figure [5] presents contours of (e 2lS ) and 4> in /J, q -T plane 
calculated by the PNJL model with scalar-type eight-quark 
interaction. Solid curves correspond to contours of (e 2ie ) = 
0, 0.4, and 0.8, while dashed curves do to contours of $ = 
0.3,0.5 and 0.7. For fi q < m^/2, the average phase fac- 
tor is finite and then calculable with LQCD in principle. For 
/i q > m w /2, (e 2ie ) = at <P < 0.5 corresponding to the 
confinement region. At <P > 0.5 corresponding to the decon- 
finement region, the factor (e ) is finite and an increasing 
function of <P, indicating that there is a positive correlation 
between (e 2ie ) and <P. Thus, LQCD is feasible only in the 
deconfinement phase, when /i q > m^/2. 
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Fig. 5: Contours of the average phase factor and the Polyakov loop 
in /iq-T plane. Contours of (e 2lf> ) — 0,0.4,0.8 and 1 are drawn 
by solid curves, while those of <P — 0.3, 0.5 and 0.7 are by dashed 
curves. 



C. Dynamical mesonic fluctuations 

The average phase factor was calculated with LQCD [6] in 
which the lattice size is 16 3 x 4 and the pion mass at vacuum 
is TO ^ QCD « 280 MeV. In the PNJL calculation, we have then 
varied the quark mass from too = 5.5 MeV to 22.5 MeV to re- 
produce tt^lQCD _ 2gQ MeV p or this value of mo, the decon- 
finement transition temperature becomes a bit higher value, 
i.e., T c = 180 MeV. 
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Fig. 6: Effects of mesonic fluctuations on the average phase factor. 
All the calculations take account of the eight-quark interaction. In 
MF+DF and MF+DF2 calculations, m n = m^ QCD and 2.8ra^ QCD , 
respectively. 

Figure [6] shows /i q dependence of the average phase fac- 
tor at T = 0.9T C ,T C and 1.25T C . LQCD result is evaluated 
from LQCD data at imaginary chemical by assuming a poly- 
nomial function. For each panel, the two solid lines delimit 



the 90% confidence level region for the extrapolation, while 
the MF+SF calculation is represented by dotted lines. The 
MF+SF calculation overestimates LQCD data largely. There- 
fore, we consider dynamical mesonic fluctuations beyond the 
MF+SF framework. 

Below T c where the system is confined, it is natural to 
think that mesonic modes dominate rather than quark modes. 
Therefore, we should take all possible channels (mesonic 
modes) in the bubble summation in the random phase approx- 
imation (RPA). If there is no pion condensate, a and 7r meson 
modes are decoupled to each other. Hence, the mesonic polar- 
ization function matrix does not have off-diagonal elements. 
Up to order 1/N C , where N c is the number of colors, the ther- 
modynamic potential is obtained by 



n = n 



MF 



Q 



DF, 



(27) 



where i?MF is the mean-field part shown in ( [Tol l and i?DF is 
the dynamical mesonic fluctuation (DF) part in the ring dia- 
gram. Hence, J7df is obtained by [44] 



■^df — — — 



d 4 q 



In det[l-G* 17(g)] 



with the effective coupling G* = diag(G* ■) for 



G* 



dM 






= G* = 



ON 
dir ' 



(28) 



(29) 



and the mesonic polarization bubbles 



n 3 k(g) 



d 4 p 



Tr[r*S(p + g)r k S(p)} (30) 



for j, k = a, 7r+, 7r_, 7To, where Tr is the trace in color, 
flavor and Dirac indexes and the four momentum integral is 
defined as J d 4 <?/(2vr) 4 = iTJ2„J d 3 q/(2?r) 3 for finite tem- 
perature. The meson vertexes 1^ depend on meson taken; pre- 
cisely, r CT = 1, r iT+ = ir+js, r 7r _ = lY_7 5] r WD = rr 3 7 5 . 

Before calculating dynamical mesonic fluctuations, we 
clarify a relation between static and dynamical mesonic fluc- 
tuations. The effective action of the NJL-type model to 1/N C 
order is derivable by using the auxiliary field method lfl7ll . 
The second derivative of the effective action with respect to 
meson fields yields the inverse meson propagator G*[l — 
G*II(q)]. Since the static fluctuations are constant, the Hes- 
sian matrix ( l25T l is obtained by setting the external momentum 
q = in the inverse propagator: 



H = G* s [l-G* s II(q = 0)). 



(31) 



Setting q = in d28l ) reduces I?df to lndet[77] and the re- 
sultant partition function turns out to be d22b . Thus, the static 
mesonic fluctuation J?sf = bidet [if] can be regarded as an 
approximation to the dynamical one 1?df- 

Now we introduce the effective meson mass to* which sat- 



isfy 



det[l - G* s n(q + fij = to*, q = 0)] = 



(32) 



with the meson chemical potentials \x a = [i Wo = 0, f/,^, = 
2fj,i and /Lt„._ = — 2/xi. Here note that physical meson masses 



to* — jij because they are calculated For /i q > m^/2, the pion condensate occurs at <P < 0.5, 



are not to* but m,j 
from q . 

Since it is difficult to calculate the dynamical mesonic fluc- 
tuations (f28b exactly, here we make the pole approximation 
that neglects the scattering phase shift. If T < T c and there 
is no pion condensation, to* is well approximated by the me- 
son mass m°, at vacuum. In this approximation, J?df can be 
obtained by a sum of four quasiparticles, a, ir + , tt~ and tt 



o. 



n 



DF 



fy 



E^> 

3 



(33) 



(34) 



where Ej — < /q 2 + to* 2 . However, corrections due to a 

and 7T° mesons are exactly cancelled out between J?i+i and 
J?i+i« . This framework is referred to as MF+DF in this paper; 
here, the static fluctuations of ^3 and 0g are taken into account 
in the Hessian matrix. 

As seen in Fig. [6] at lower T such as T = 0.9T C and T c , the 
MF+DF calculation (dashed line) almost reproduces LQCD 
data. Above T c such as T = 1.25T C , the MF+DF calcula- 
tion underestimates LQCD data. For T > T c , in general, the 
pion mass becomes larger than m° H37I1 . If m ff = 2.8to° 
is taken, the calculation (dot-dashed line) is consistent with 
LQCD data; this calculation is denoted by MF+DF2 in Fig. [6] 



IV. SUMMARY 

We have calculated the average phase factor of the QCD 
determinant at finite quark chemical potential /i q , using the 
two-flavor version of the PNJL model with the scalar-type 
eight quark interaction, since the model is successful in re- 
producing LQCD data not only on the 1+1 system in the limit 
of no /i q 1 35] but also on the 1+1* system with finite isospin 
chemical potential /ii [38]. In the present model, there exists a 
critical endpoint (CEP) in the 1+1 system with finite /i q . The 
CEP lies inside the (e 2lB ) = region. This implies that the 
location cannot be determined by LQCD solely. 



so that (e 



2/»\ 



= there. At <P > 0.5, the factor (e 



2itt\ 



IS 



finite and an increasing function of <P. Thus, there exists a 
positive correlation between (e 2,10 ) and <P. Therefore, LQCD 
is feasible only in the deconfinement phase with large 4>, when 

/iq > TO^/2. 



The eight-quark interaction makes the 



= region 



shrink, while it shifts the CEP toward higher T and smaller 
/i q . As a consequence of these properties, a relative distance 
of the CEP to the boundary of the (e 2l °) = region becomes 
smaller. If more accurate LQCD data becomes available in 
future outside the (e 2l °) = region, we can predict a location 
of the CEP with the PNJL model the parameters of which are 
fitted to the data. The accuracy of the model prediction seems 
to be proportional to the distance. In this sense, the fact that 
the relative distance is small is important. 

For /x q < TO,r/2 where no pion condensate takes place, 
the PNJL calculation with static fluctuations cannot reproduce 
LQCD data @] at both T < T c and T > T c . This problem is 
solved partly by treating dynamical mesonic fluctuations with 
the pole approximation. The PNJL model with the dynamical 
mesonic fluctuations reproduces LQCD data at T < T c , For 
T > T c , however, the calculation cannot reproduce the data. 
The first possible reason is that the meson mass at T > T c is 
different from the value at vacuum. The second possible rea- 
son is that the pole approximation is not good, because meson 
at T > T c is generally in a resonance state. The third possible 
reason is effects of physical states with nonzero baryon num- 
bers. It is reported in Ref. [6] that the physical states tend to 
enhance the average phase factor, although the PNJL model 
does not include such effects. This problem should be solved 
in future in order to construct a reliable effective model that 
makes it possible to predict a location of CEP and the phase 
diagram at finite /u q . 
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